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^ ; Abstract 

"q ■ We investigated the cosmology in a higher-curvature gravity where the 

dimensionality of spacetime gives rise to only quantitative difference, 
^ . contrary to Einstein gravity. We found exponential type solutions for 

flat isotropic and homogeneous vacuum universe for the case in which 
^ \ the higher- curvature term in the Lagrangian density is quadratic in the 

scalar curvature, £i? 2 . The solutions are classified according to the sign 
of the cosmological constant, A, and the magnitude of A£. For these 



in 



solutions 3-dimensional space has a specific feature in that the solutions 



o 



are independent of the higher curvature term. For the universe filled with 
perfect fluid, numerical solutions are investigated for various values of the 
parameter £. Evolutions of the universes in different dimensionality of 
^ | spacetime are compared. 

SqI PACS numbers: 04.50.+h, 98.80.-k 



1 Introduction 



Einstein gravity describes the observed universe fairly well. It predicts, however, the presence 
of initial singularity in our universe This is usually interpreted to mean the limit 
of applicability of the Einstein gravity. Candidates for the alternative theory come from 
both classical and quantum considerations. In classical regime, there exists a possibility of 
avoiding the initial singularity in higher-curvature gravity(HCG) theories. HCG seems to 
be natural when the gravity is very strong and the curvature is large, since linearity in R of 
the Lagrangian would be too restrictive. Numerical investigations of the initial singularity 
in HCG have already been carried out [3] and singularity free solutions have been found, 
although the solutions are limited only to a very early stage of the universe. HCG is also 
suggested from the consistency of quantum theory, e.g. quantum field theory in curved 
spacetimes or string perturbation theory. In these cases tidal effects are important and 
equivalence principle would not work effectively. 

Thus HCG seems to work in both classical and quantum theory. From the cosmological 
point of view, it may give a better description of the initial stage of the universe than Einstein 
one and is expected to give a similar description of the later universe as the Einstein one 
since the curvature effects would then be small. However, the whole evolution of the universe 
in HCG has not been investigated in detail so far and the description of the present state of 
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the universe is not satisfactory jlj. 

There have been many investigations on cosmology in superstring motivated HCG theo- 
ries. They can be classified into two types. In one of them, the results of the perturbation 
theory are used and the Gauss-Bonnet combination is adopted. In another, a suitable 
hypothes is adopted to construct HCG Lagrangian density leading to produce non-singular 
cosmological solutions, instead of using the results of the perturbation theory. [SJ[71IBJIII| Both 
types of models yield interesting cosmological solutions, although they are not satisfactory, 
since the solutions do not describe the present state of the universe. 

Theoretically the Gauss-Bonnet combination has a nicely simple property. However it is 
so simple that new generalized coordinates, coming from the time derivatives of the metric 
and inherent to HCG, do not appear in the canonical formalism at least for the Robertson- 
Walker spacetime. Therefore it would be worthwhile to investigate other types of HCG. 
One of the simple models is described by a function of the scalar curvature. This type of 
models has long been investigated and includes many problems peculiar to HCG such as 
the definition of the new generalized coordinate, complicated constraint structure PUI IT2^ IT3] 
and conformal equivalence Thus it is interesting to see whether this type of models are 
applicable to describe the realistic evolution of the universe. 

In this paper, we investigate the evolution of the universe in HCG, described by a function 
of the scalar curvature, using analytical and numerical methods. In numerical analyses, 
we adopt the recent observational data as the initial data which leads to the accelerated 
expansion, instead of assuming suitable initial data at t = 0. One of the main differences of 
HCG from Einstein one is that the dimension of the spacetime gives rise to only quantitative 
differences, e.g. while in Einstein gravity there is no dynamical degrees of freedom in (2+1)- 
dimensional spacetime, there remains such freedom in HCG ^21 QHI- Therefore it is easy 
to compare spacetimes differing in dimensionality [T%] . In other words it would be easy to 
transfer results obtained for lower dimensional spacetime to higher dimensional one. We 
investigated the universe for the case of vacuum universe and the universe filled with perfect 
fluid for various values of the coefficient £ of the i? 2 -term. It is found that (3+l)-dimensional 
spacetime has somewhat peculiar property in the former case. 

In section 2 we present the basic equations for general spacetimes as well as for homoge- 
neous and isotropic universe in the case of R 2 gravity. In section 3 analytical and numerical 
solutions are investigated. We obtain exponentially evolving solutions for a particular model 
containing a term quadratic in scalar curvature. The results of numerical investigations for 
the vacuum case and the universe filled with perfect fluid are presented for various values of 
£. The solutions exhibit largely different behavior according to the sign of £. Section 4 is 
devoted to the summary and discussions. 

2 Formulation 

2.1 Basic equations 

We consider a pure gravity type model in (d + l)-dimensional spacetime described by a 
Lagrangian density 



where / is a differentiable function of the (d + l)-dimensional scalar curvature R and a hat 
denotes a quantity defined in (d + l)-dimensional spacetime. Field equations derived from 
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(2.1) is known to be written as 

f'R, v - l -g, v f + g^Uf - V M V,/' = (2.2) 

where R pv is the Ricci tensor constructed from the (d + l)-dimensional metric g^, V M is the 
covariant derivative with respect to g^ u , □ = g xp V\V p and /' = df /dR. 

Equations (2.2) depend on the dimension of spacetime only implicitly. The dependence 
becomes explicit in the canonical formalism. The system described by the Lagrangian density 
(2.1) is highly constrained, so we follow the method of Buchbinder and Lyakhovich (BL) ^0} 
12j for treating the higher-derivative Lagrangian. The Hamiltonian density 7i* is obtained 
by the Legendre transformation of the modified Lagrangian density C* G which takes into 
account the constraints and the definitions of the generalized coordinates coming from the 
time derivatives of the metric by introducing the Lagrange multipliers. 7i* has the following 
form: 

H* G = NH G + N k H k + divergent terms (2.3) 

where N is the lapse function and N h is the shift vector. Tic — and Tik = are the Hamil- 
tonian and momentum constraints, respectively. Choosing the coordinate system where the 
lapse function N is equal to 1 and the shift vector N h vanishes, explicit form for Tic — 7~L*g 
is given as follows [T2"]: 

Kg = 2U- 1 (jPpy - ^p 2 ) + ^ P Q + hlR - Vhf(R) - ^HQ 2 + An - l -im (2.4) 

where hij is the metric of the <i-dimensional space, p lJ the momentum canonically conjugate 
to h^, h = det h^, 1Z the scalar c urvature formed from hy. p = h^p 11 , Q = h^Q^j where 
Qij is the extrinsic curvature which is taken as new generalized coordinate replacing the 
time-derivative of hij. IT is the momentum canonically conjugate to Q. The scalar curvature 
of the (d + l)-dimensional spacetime, R, should be expressed in terms of the canonical 
variables through a relation II = 2\/hf'(R). It is seen that the dimension d appears only in 
the coefficients of some of the terms in such a way that there is no special value of d. 

2.2 Homogeneous, isotropic and flat spacetimes 

In order to solve the field equations explicitly, we have to specify the model. Here we 
specialize the function f(R) to a quadratic function: 

f(R) = -2A + R + £R 2 (2.5) 

Then, in terms of canonical variables, R is expressed as 

R= ^(Vn/v^-l) (2.6) 

where k = 8irG. Furthermore we specialize the spacetime to the flat Robertson- Walker 
spacetime: 

ds 2 = -dt 2 + a(t) 2 h ij dx i dx j (2.7) 

where is the metric of flat (i-dimensional static space. Then H,g in (2.4) reduces to the 
following form: 

H G = ^a- 2 U 2 - I ( 3 Q 2 + i) H + \ Va aQ + ~ (a + 1) a 2 (2.8) 
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(2.9) 



where p a is the momentum canonically conjugate to the scale factor a. 

For vacuum spacetime the canonical equations of motion are derived from (2.8): 

A => p-w a " (J+I,n2 ->^( A+ 4) a "- 

q = ^a-"n - 1 (^±iq 2 + 1), fi = - d {(d + i)qn - a Pa ) 

From these equations, we have an equation for a 1 : 

(d - 1)- + J(d - l)(d - 2) (-V - A + ifCT = 

where HCT represents the contribution of the higher curvature terms and is expressed as 

,( 4 K /Ax /n .(3) N ^ N 2 



(2.10) 



ifCT 



d£ 



*(v) + «--^©(v) + w-^e)' 



(2.11) 



+2(3d 2 -20d + 21)(-) (-) +-(d-l)(d-4)(d-9)(-) 



3 Solutions 

3.1 Analytical solutions 

It is easily seen that (2.10) allows a particular solution of the form a(t) = a e xt with constant 
a . Putting this form into (2.10), we have 

d(d - 1)A 2 - 2A + £d 2 (d + l)(d - 3)A 4 = (3.1a) 

so that 



A 2 = -d(d - 1) ± ^(d - l) 2 + 8A£d 2 (d + l)(d - 3) /2£d 2 (d + l)(d-3) (3.1b) 

Thus the solutions are classified according to the value of 

?? = 8A£(d+l)(d-3)/(d-l) 2 , 
and the sign of A. The solutions exist only for i] > — 1. 



(3.2) 



I. a(t) = a Q exp 



II. a(t) = ao exp 



III. a(t) = a exp 



a(t) = do exp 



±1 



4A 



d(d- 1) 



t 



for 7] = — 1 



and A > 



d(d — l)r] 



for - 1 < 7] < and A > 



± 



\Z*r^(- 1 + ^) ( 



for 77 > 



for 77 > 



and A > 

and A < 
(3.3) 



1 Of course, (2.10) can be derived from the space components of (2.2). 
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This classification can be made more directly (although not compact), e.g. 77 = — 1 means 
Af = -(d - l) 2 /8(d + l)(d - 3) or r] > means A£(d - 3) > 0. It should be noted that, 
for d = 3, (3.1a) shows that there is no contribution from the higher curvature term, so the 
solution is the same as in Einstein gravity. This is also seen if the trace of (2.2) is taken 
and (2.5) is used for f(R). Correspondingly, for d — 3, and consequently 77 = 0, only half 
of the type II and type III solutions have i] —>■ limit. We note that, for A = 0, nontrivial 
solutions exist only for d = 2 if £ > 0. It would be expected that these exponential type 
vacuum solutions are relevant only to the early stage of the universe. Thus, in order to 
examine whether the initial singularity of our universe can be avoided or not, other types of 
solutions are necessary. They will be investigated numerically in the next subsection. 

In later stage of the universe, we assume the universe to be filled with perfect fluid, when 
the field equations are obtained by replacing the rhs of (2.2) by the energy-momentum tensor 
of the perfect fluid. Equation for the scale factor a is expressed as 



U - l)_ + -(d - l)(d - 2) (-Y - A + HCT = k P . (3.4) 
a 2 \aJ 

The field equations are supplemented by the equation of state 

P = IP ( 3 - 5 ) 

where p is the pressure, p the total energy density and 7 a constant. Energy-momentum 
conservation leads to a relation: 

/ a \ ^(7+1) 

P = Po - ■ 3.6 

Va / 

Analytical solutions are available only for radiation dominated era, 7 = 1/0?, which reads 

a(t) oc t 2/l ~ d+1) . (3.7) 

Also in this case, d = 3 is exceptional in that the solution (3.7) is the same as in Einstein 
gravity. Other solutions are investigated numerically in the next section. 



3.2 Numerical analysis 

In this section we present the results of the numerical analyses of the equation (2.10), for the 
vacuum case, and (3.4) with (3.5), describing the universe filled with perfect fluid. These 
equations were solved numerically by the 4th order Runge-Kutta method. 

We start from the universe today and investigate the evolution of the universe both to 
the future and to the past. The present values of the cosmological parameters are taken from 
the WMAP results [15J. The energy densities of matter pmo and dark energy (cosmological 
constant) pdeo are taken to satisfy 

(i) pmo ■ Pdeo = 27 : 73, and (ii) Q M o + ^deo = 1 (flat universe). (3.8) 

The latter means, of course, the total energy density is critical. The Hubble constant is 
taken to be 

H = 72kms- 1 Mpc" 1 . (3.9) 

(3) 

The value of do is determined approximately from the Einstein equation and a from the 
Hamiltonian constraint. It is noted that these values of parameters lead to the accelerating 
universe as is shown in the following figures. 
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Figure 1 shows the predicted evolutions of the universe filled with perfect fluid for d = 3 
according to Einstein gravity and HCG for various values of the parameter £. For positive £, 
the evolutions predicted by HCG are almost the same as that by Einstein gravity as expected. 
On the other hand, for negative £, predictions of HCG are very different from that of Einstein 
gravity and depend heavily on the values of £. These results remain unchanged for the value 
of Hq other than the one in (3.9). 2 . 

Figure 2 shows the evolutions of the scalar curvature for the cases as in Figure 1. The 
curvature becomes very large at early stage as expected and corresponds to the behavior of 
the scale factor. The results shown in Figures 1 and 2 seem to rule out the cases of negative 
£ since the calculated age of the universe is so short. 

Figure 3 shows the evolution of the scale factor for the dimensions of the space, d = 2, 3, 4 
in Einstein gravity. As to the initial conditions for d = 2 and d = 4, we tentatively adopted 
the same ones as for d = 3. Figure 4 shows the corresponding evolution of the scalar 
curvature. 

Figures 5 and 6 show the same evolution as in Figures 3 and 4 in HCG for £ = 10. 
The difference in the equation of motion is only quantitative as noted above, which is in 
accordance with our numerical results. We note that our numerical solutions satisfy the 
Hamiltonian constraint fairly well. 




x 10° yr 



Fig.l Evolutions of the (3+l)-dimensional universe filled with perfect fluid 
in Einstein gravity and HCG for various values of £ 
2 The older values of these parameters |lfil IT7I I18I fTTH I2()j lead to almost the same results. 
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Fig. 2 Evolutions of scalar curvature R in (3+l)-dimcnsional spacetime 





Fig. 3 Evolutions of the universe in Einstein grav- 
ity for dimensions d — 2,3 and 4. 



Fig. 4 Evolutions of scalar curvature R in Einstein 
gravity for dimensions d — 2, 3 and 4. 
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x 1 s yr x10 8 yr 

Fig.5 Evolution of the universe in HCG with £ = Fig.6 Evolutions of scalar curvature R in HCG 

10, for dimensions d = 2, 3 and 4 for £ = 10 for dimensions d = 2, 3 and 4 



4 Summary and discussion 

We investigated the evolution of the flat, RW universe in a higher curvature gravity (HCG) 
theory described by a Lagrangian depending on the scalar curvature. Dimensionality of 
spacetime affects the evolution only quantitatively. Explicit solutions are obtained when 
the higher curvature term is quadratic in the scalar curvature, £i? 2 . Solutions exhibiting 
exponential expansion are obtained for the vacuum universe. Exceptional case is that of 
3-dimensional space, in which the solution is the same as in Einstein gravity, i.e. the higher 
curvature term gives rise to no effect. These solutions may be interpreted to mean that 
other type of solution is necessary to avoid the initial singularity since the universe has not 
expanded exponentially throughout its evolution. 

Numerical investigations were made by taking the recent observational data as initial 
conditions. The solutions were compared for various values of the parameter £ and Einstein 
gravity. For positive £, solutions in HCG deviate only slightly from that in Einstein gravity. 
For negative £, solutions depend heavily on the value of £ and indicate too short an age 
of the universe. However, numerical results do not show indications of avoiding the initial 
singularity in the range where our numerical calculations are valid. It may be that a single 
set of equations is not sufficient to describe the whole history of the universe and different 
set of equations are required for the earliest and later stages of the universe analogous to 
the usual inflationary scenario. In other words, it may be required that numerical solutions 
should be connected to solutions, e.g. those obtained in the previous section at some early 
time. 

Concerning the initial singularity, it is known that for spacetime with physically rea- 
sonable properties, occurrence of the singularity is inevitable if the Ricci tensor satisfies 
RavC^C > f° r an Y timelike vector In HCG with £i? 2 term, we have 



1 + 2ZR 



1 



d-l 



T 



d-l 



A 



(4.1) 



(d-l)(l + 2£i2)L 



R 2 + 2DR - 2(d - 1)CCV^ U R 



The second term on the rhs is due to the higher curvature effect. It is difficult to determine 
the sign of the rhs generally, even if the strong energy condition is satisfied. It is seen that 
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the results may change according the sign of £. Restricting to the isotropic and flat universe 
filled with pressureless perfect fluid, we have, in terms of the canonical variables 



V 




1 



2 



A 



d-l 



p 



d- 1 



1 



k 2 T1 Vh 



- 2 



□n 



2(d-i)e 



K 2 n 



(d- i)n 



+ 



n 



Even in this case the situation is hardly improved so that numerical analysis seems to be 
necessary. Our numerical results show that R takes large positive value for £ > 0. This 
indicates the possibility that the second term on the rhs of (4.1) dominates the first term 
and the numerical solutions could be connected to some non-singular solution. 

We also compared the evolution of the universe for the dimensions of space d = 2, 3 and 4. 
Intuitively the case of lower dimensional spacetime would seem to be simpler |14j . However 
we found that the dimension of the spacetime is almost irrelevant in HCG. 

Finally we comment on the conformal transformations in HCG^T]. In this work we iden- 
tified the metric g^ v in (2.1) as physical and applied the initial conditions to them. If the 
transformed metric is identified to be physical, the initial conditions are changed, although 
the equations of motion are only transformed. The transformation is suggested by the form 
of Tic i n (2.4) which is so complex. The form can be brought to a simpler form. Investiga- 
tion of the transformed equations would be interesting since new types of solutions could be 
obtained. 
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